Theory of Frozen Waves ^ 



M. Zamboni-Rached, 
DMO-FEEC, State University at Campinas, Campinas, SP, Brazil. 

Erasmo Recami 

Facolta di Ingegneria, Universita statale di Bergamo, Dalmine (BG), Italy; 
and INFN — Sezione di Milano, Milan, Italy. 

and 

H. E. Hernandez-Figueroa 
DMO-FEEC, State University at Campinas, Campinas, SP, Brazil. 

Abstract — In this work, starting by suitable superpositions of equal-frequency Bessel 
beams, we develop a theoretical and experimental methodology to obtain localized sta- 
tionary wave fields, with high transverse localization, whose longitudinal intensity pattern 
can approximately assume any desired shape within a chosen interval < z < L of the 
propagation axis z. Their intensity envelope remains static, i.e. with velocity v = 0; so 
that we have named "Frozen Waves" (FW) these new solutions to the wave equations 
(and, in particular, to the Maxwell equations). Inside the envelope of a FW only the 
carrier wave does propagate: And the longitudinal shape, within the interval < z < L, 
can be chosen in such a way that no nonnegligible field exists outside the pre-determined 
region (consisting, e.g., in one or more high intensity peaks). Our solutions are noticeable 
also for the different and interesting applications they can have, especially in electromag- 
netism and acoustics, such as optical tweezers, atom guides, optical or acoustic bistouries, 
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1 Introduction 



Over many years the theory of localized waves (LW), or nondiffracting waves, and in 
particular of the so-called X-shaped waves, has been developed, generalized, and exper- 
imentally verified in many fields, such as optics, microwaves and acoustics lj. These 
new solutions to the wave equations (and, in particular, to the Maxwell equations) have 
the noticeable characteristic of resisting the diffraction effects for long distances, i.e., of 
possessing a large depth of field. 

Such waves can be divided into two classes: the localized beams and the localized 
pulses. With regard to the beams, the most popular is the Bessel beam. 

Much work was made about the properties and applications of a single Bessel beam, 
while some work has been done in connection with Bessel beam superpositions performed 
by summing or integrating over their frequency (by producing, e.g., the well-known "X- 
shaped pulses" and/or their velocity (for example, it has been studied the space-time 
focusing of different-speed X-shaped pulses). By contrast, only a few papers have been 
addressed to the properties and applications of superpositions of Bessel beams with the 
same frequency, but with different longitudinal wave numbers. The few works existing 
on this subject have shown some surprising possibilities associated with this particular 
type of superpositions, mainly the possibility of controlling the transverse shape of the 
resulting beamj2El- 

The other important point, i.e., that of controlling the longitudinal shape, has been 
even more rarely addressed, and the relevant papers have been so far confined to numerical 
optimization processes jH Ej, aimed at finding out an appropriate computer-generated 
hologram. 
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In this work we develop a very simple method** , having recourse to superpositions 
of forward propagating and equal-frequency Bessel beams only, that allows controlling 
the beam-intensity longitudinal shape within a chosen interval < z < L, where z is the 
propagation axis and L can be much greater than the wavelength A of the monochromatic 
light (or sound) which is being used. Inside such a space interval, indeed, we succeed in 
constructing a stationary envelope whose longitudinal intensity pattern can approximately 
assume any desired shape, including, for instance, one or more high-intensity peaks (with 
distances between them much larger than A); and which results — in addition — to be 
naturally endowed also with a good transverse localization.*** 

This intensity envelope remains static, i.e., has velocity v = 0; and because of this in 
a previous paper we have called "Frozen Waves" (FW) these new solutions to the wave 
equations (and, in particular, to the Maxwell equations). Inside the envelope of a FW 
only the carrier wave does propagate: And the longitudinal shape, within the interval 
< z < L, can be chosen in such a way that no nonnegligible field exists outside the 
pre-determined high-intensity region. 

We also suggest a simple apparatus capable of generating the mentioned stationary 
fields. 

Static wave solutions like these are noticeable also for the different and interesting 
applications they can have, especially in electromagnetism and acoustics, such as optical 
tweezers, atom guides, optical or acoustic bistouries, optical micro-lithography, electro- 
magnetic or ultrasound high-intensity fields for various important medical purposes, and 
so on.** 



** Patent pending. 

*** When we get a complete control on the longitudinal shape, we cannot have — however — a total 
control also on the transverse localization, since our stationary fields are of course constrained to obey 
the wave equation. 
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2 The mathematical methodology:** Stationary 
wavefields with arbitrary longitudinal shape, ob- 
tained by superposing equal-frequency Bessel 
beams 

Let us start from the well-known axis-symmetric zeroth order Bessel beam solution to the 
wave equation: 

r(>(j>,z,t) = J (k p p)e^e~^ (1) 

with 

k>^~P\ (2) 

where u, k p and (3 are the angular frequency, the transverse and the longitudinal wave 
numbers, respectively. We also impose the conditions 

u/P > and k 2 p > (3) 

(which imply uj / (3 > c) to ensure forward propagation only (with no evanescent waves), 
as well as a physical behavior of the Bessel function Jo- 

Now, let us make a superposition of 2N + 1 Bessel beams with the same frequency 
ojq, but with different (and still unknown) longitudinal wave numbers /3 n : 

N 

*{p,z,t) = e"'^* £ A n J Q (k pn p)e ll3 " z , (4) 

n=-N 

where n are integer numbers and A n are constant coefficients. For each n, the parameters 
ljq, k pn and (3 n must satisfy Eq.(J2J), and, because of conditions ©, when considering 
ujq > 0, we must have 
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Let us now suppose that we wish \^f(p, z, t)\ 2 of Eq.Q to assume on the axis p = 
the pattern represented by a function |F(z)| 2 , inside the chosen interval < z < L. In 
this case, the function F(z) can be expanded, as usual, in a Fourier* series: 

oo 

F{z) = £ B m e^ mz , 

m=—oo 

where 

B m = \ f L F(z)e-^ mz dz . 

More precisely, our goal is finding out, now, the values of the longitudinal wave numbers 
f3 n and of the coefficients A n , of Eq.( | % | ). in order to reproduce approximately, within the 
said interval < z < L (for p = 0), the predetermined longitudinal intensity-pattern 
|.F(,2)| 2 . Namely, we want to have 

2 

w \F(z)\ 2 with < z < L . (6) 

Looking at Eq.(JBJ), one might be tempted to take (3 n = 27m/ L, thus obtaining a 

truncated Fourier series, expected to represent approximately the desired pattern F(z). 

*Such a choice of the longitudinal intensity pattern does imply an interesting freedom, since we can 
consider more in general any expansion X)m=-oo ^ m ex P* ~r~ mz = F(z) expi0(z), quantity 4>{z) being 
an arbitrary function of the coordinate z. 



N 
n=-N 
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Superpositions of Bessel beams with f} n = 2im/ L have been actually used in some works 
to obtain a large set of transverse amplitude profiles [21 E]. However, for our purposes, this 
choice is not appropriate, due to two principal reasons: 1) It yields negative values for 
f3 n (when n < 0), which implies backwards propagating components (since uj > 0); 2) 
In the cases when L » Ao, which are of our interest here, the main terms of the series 
would correspond to very small values of /3 n , which results in a very short field-depth 
of the corresponding Bessel beams (when generated by finite apertures), preventing the 
creation of the desired envelopes far form the source. 

Therefore, we need to make a better choice for the values of f3 n , which allows forward 
propagation components only, and a good depth of field. This problem can be solved by 
putting 



where Q > is a value to be chosen (as we shall see) according to the given experimental 
situation, and the desired degree of transverse field localization. Due to Eq.(JSJ), we get 



Inequality (jHJ), can be used to determine the maximum value of n, that we call iV max , 
once Q, L and uiq have been chosen. 

As a consequence, for getting a longitudinal intensity pattern approximately equal to 
the desired one, |F(2;)| 2 , in the interval < z < L, Eq.(@J) should be rewritten as 



Pn = Q + T" n ' 



(7) 



2 71 UOn 

< Q± — N < — . 

L c 



(8) 



N 



*(p = 0,z,t) 



e -iu t e iQz ^ A n e 



n z 



(9) 



n=-N 



with 




(10) 
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Obviously, one obtains only an approximation to the desired longitudinal pattern, because 
the trigonometric series © is necessarily truncated (N < iV max ). Its total number of 
terms, let us repeat, will be fixed once the values of Q, L and ujq are chosen. 
When p ^ 0, the wavefield *&(p, z, t) becomes 

2V 

*(p,z,t) = e~^ t e iQz A n J Q (k pnP )e^ nz , (11) 

n=~N 

with 

Kn = ^ -(Q + 2 -^) 2 . (12) 

The coefficients A n will yield the amplitudes and the relative phases of each Bessel 
beam in the superposition. 

Because we are adding together zero-order Bessel functions, we can expect a high field 
concentration around p = 0. Moreover, due to the known non-diffractive behavior of the 
Bessel beams, we expect that the resulting wavefield will preserve its transverse pattern 
in the entire interval < z < L. 

The methodology developed here deals with the longitudinal intensity pattern control. 
Obviously, we cannot get a total 3D control, due the fact that the field must obey the wave 
equation. However, we can use two ways to have some control over the transverse behavior 
too. The first is through the parameter Q of Eq.(7). Actually, we have some freedom 
in the choice of this parameter, and FWs representing the same longitudinal intensity 
pattern can possess different values of Q. The important point is that, in superposition 
using a smaller value of Q makes the Bessel beams possess a higher transverse 
concentration (because, on decreasing the value of Q, one increases the value of the Bessel 
beams transverse numbers), and this will reflect in the resulting field, which will present 
a narrower central transverse spot. We will exemplify this fact in the next Section. The 
second way to control the transverse intensity pattern is using higher order Bessel beams, 
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but we shall show this in Section 5. 



3 Some examples 

In this Section we shall present a few examples of our methodology. 
First example: 

Let us suppose that we want an optical wavefield with A = 0.632 /im, that is, with 
ujq = 2.98 x 10 15 Hz, whose longitudinal pattern (along its z-axis) in the range < z < L 
is given by the function 



F(z) 



(z-h){z-l 2 ) 
(h-h) 2 

1 

(z - k)(z - Z 6 ) 
(k ~ h) 2 





for l\ < z < I2 
for I3 < z < U 
for l 5 < z < Iq 
elsewhere , 



(13) 



where l\ = L/5 — Az±2 and I2 = L/5 + Azu with Az\2 = L/50; while Z3 = L/2 — Az^ and 
U = L/2 + A^ 3 4 with Az 34 = L/10; and, at last, I5 = 4L/5 — Az 56 and Iq = 4L/5 + Az 56 
with A^ 56 = L/50. In other words, the desired longitudinal shape, in the range < z < L, 
is a parabolic function for h < z < l 2 , a. unitary step function for / 3 < z < / 4 , and again a 
parabola in the interval h < z < Iq, it being zero elsewhere (within the interval < z < L, 
as we said). In this example, let us put L = 0.2 m. 

We can then easily calculate the coefficients A n , which appear in the superposition 
(JTTJ), by inserting Eq.(JT3J) into Eq.tJTHJ). Let us choose, for instance, Q = 0.999 ujq/c: This 
choice allows the maximum value iV max = 316 for n, as one can infer from Eq.fjHj). Let 
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us emphasize that one is not compelled to use just iV = 316, but can adopt for N any 
values smaller than it; more generally, any value smaller than that calculated via Eq.(JBJ). 
Of course, on using the maximum value allowed for N, one gets a better result. 

In the present case, let us adopt the value N = 30. In Fig. 1(a) we compare the 
intensity of the desired longitudinal function F(z) with that of the Frozen Wave, $?(p — 
0, z, t), obtained from Eq.Q by adopting the mentioned value N = 30. 



hi | | 
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0.6 ■ ! i 

0.4 ■ 
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z[m] 

(a) 

Figure 1: (a) Comparison between the intensity of the desired longitudinal function F(z) 
and that of our Frozen Wave (FW), \I/(p = 0, z,t), obtained from Eq. ([§]). The solid line 
represents the function F(z), and the dotted one our FW. (b) 3D-plot of the field- intensity 
of the FW chosen in this case by us. 

One can verify that a good agreement between the desired longitudinal behavior and our 
approximate FW is already got with iV = 30. The use of higher values for N can only 
improve the approximation. Figure 1(b) shows the 3D- intensity of our FW, given by 
Eq. (lll|) . One can observe that this field possesses the desired longitudinal pattern, while 
being endowed with a good transverse localization. 

Second example: 
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Let us now suppose we want an optical wavefield with A = 0.632 /mi (uq = 2.98 1015 Hz), 
whose longitudinal pattern (on its axis) in the range < z < L consists in a pair of 
parabolas, for l\ < z < l 2 and h < z < U, the intensity of the second parabola being 
twice as much as that of the first one. Outside the intervals l\ < z < l 2 U h < z < U, 
the desired field has a null intensity. Summarizing, we want : 

-4 {Z - h)[Z - 2 h) for h<z<h 



F(z) 



-4^2 {z - h){z - h) for l 3 <z< h ™ 
[k - h) 



elsewhere , 

where h = 3L/10-Az 12 and l 2 = 3L/10+Az 12 with Az 12 = L/70; while l 3 = 7L/10-Az u 
and I4 = 7L/10 + A^ 34 with Az 34: = L/70. In this example we choose L = 0.02 m. 

Again, we can calculate the coefficients A n by inserting Eq. (fT^|) into Eq.(|l()|l. and use 
them in our superposition @. In this case, we chose Q = 0.995 uj /c: This choice allows 
a maximum value of n given by iV max = 158 (one can see this by exploiting Eq.(jHJ). But 
for simplicity we adopt once more N = 35, hoping that Eq.(|lip will yield a good enough 
approximation of the desired function. 

We compare in Fig.2(a) the intensity of the desired longitudinal function F(z) with 
that of our FW, \&(p = 0, z, t), obtained from Eq.(JH} by using iV = 35: We can verify a 
good agreement between the desired longitudinal behaviour and our FW. Obviously we 
can improve the approximation by using larger values of N. 

In Fig. 2(b) we show the 3D field intensity of our FW, forwarded by Eq. ljllj) . We can 
see that this field has a good transverse localization and possesses the desired longitudinal 
pattern. 
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Figure 2: (a) Comparison between the intensity of the desired longitudinal function F(z), 
given by Eq.([14|). and that of our FW, ^(p = 0, z, t), obtained from Eq.(J2J). The solid 
line represents the function F(z), and the dotted one our FW.(b)3D plot of the field 
intensity of the FW chosen by us in this new case. 

Third example (controlling the transverse shape too): 
We want to take advantage of this new example for addressing an important question: 
We can expect that, for a desired longitudinal pattern of the field intensity, by choosing 
smaller values of the parameter Q one will get FWs with narrower transverse width [for 
the same number of terms in the series entering Eq. ljllj) ]. because of the fact that the 
Bessel beams in Eq. ljllj) will possess larger transverse wave numbers, and, consequently, 
higher transverse concentrations. We can verify this expectation by considering, for 
instance, inside the usual range < z < L, the longitudinal pattern represented by the 
function 

-4 V K ;- fa) ^ h <z<l 2 

{h-li) ? (15) 

elsewhere 

with l\ — L/2 — Az and 1% = L/2 + Az. Such a function has a parabolic shape, with its 
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peak centered at L/2 and with longitudinal width 2A^/v / 2- By adopting Ao = 0.632 fim 
(that is, uj = 2.98 x 10 15 Hz), let us use the superposition (fTTj) with two different values 
of Q: We shall obtain two different FWs that, in spite of having the same longitudinal 
intensity pattern, will possess different transverse localizations. Namely, let us consider 
L = 0.06m and Az = L/100, and the two values Q = 0.999 a;o/c and Q = 0.995 a;o/c. 
In both cases the coefficients A n will be the same, calculated from Eq.([10p. on using this 
time the value N = 45 in superposition fTTj) . The results are shown in Figs. 3 (a) and 3(b). 
Both FWs have the same longitudinal intensity pattern, but the one with the smaller Q 
is endowed with a narrower transverse width. 



p I 2 




Figure 3: (a) The Frozen Wave with Q = 0.999cjo/ c an d = 45, approximately repro- 
ducing the chosen longitudinal pattern represented by Eq.(|15|). (b) A different Frozen 
wave, now with Q = 0.995co>o/c (but still with iV = 45) forwarding the same longitudinal 
pattern. We can observe that in this case (with a lower value for Q) a higher transverse 
localization is obtained. 

In Section 5 we shall show that a better control of the transverse shape can be obtained 
by using higher order Bessel beams in superposition (JTTJ). 
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4 Spatial resolution and Residual intensity 

In connection with a FW of a given frequency, it is of practical (and theoretical) 
interest to investigate its Spatial resolution, its Residual intensity, the Size of the source 
necessary to generate it, as well as the minimum distance from the source needed to get 
such a FW. 

Let us first comment that, in lossless media, the theory of FWs can furnish results 
similar to the free-space ones. This happens because FWs are suitable superpositions 
of Bessel beams with the same frequency, so that there is no problem with the material 
dispersion. 

Here, we deal with lossless media only. 

Let us address the question of the longitudinal and transverse spatial resolution for 
the FWs. 

In connection with the longitudinal case, once we choose a desired longitudinal inten- 
sity field configuration, |F(z)| 2 , given, for example, by a single peak (or a few peaks) with 
a certain longitudinal width Az, we wish to investigate whether it is possible to obtain 
such a spatial resolution, and what are the relevant parameters for getting good results. 

As one can expect, this question is directly related to the number 2N + 1 of terms in 
superposition (fTTj) : More specifically, in superposition 0. 

Once the values of the frequency u , and the parameters L and Q are chosen, the 
best approximation that we can get for a given longitudinal intensity-field configuration, 
|F(^)| 2 , is obtained by using Eq.(JHJ) with the maximum number of terms 2iV max -|- 1, where 
N max is calculated from is calculated through inequality (JEJ). 

As we have seen in the previous Sections, it is not always necessary to use N = iV max , 
and frequently a smaller value of N can provide us with good results. But even in this 
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cases, when a value N < iV max is quite sufficient to furnish the desired spatial resolution, 
nevertheless it can be desirable to increase the value of N for lowering the longitudinal 
residual intensities, as we are going to see. 

In the cases in which not even the value N = iV max yields a good result, we have to 
adopt a smaller value for the parameter Q so to increase, in this way, the value of N max 
itself. For quantifying mathematically the precision of our approximation, one may have 
recourse to the mean square deviation, D, 

N 



D = [ \F(z)\ 2 dz-Lj2\A n \ 2 , 



where A n , the coefficients of superposition (jHJ), are given by Eq. (finj) . 

The case of the transverse spatial resolution cannot be tackled in such a detail, since, 
as we know, one cannot have a complete three-dimensional control of the field. In the 
previous Section we have seen that we can get, however, some control on the transverse 
spot size through the parameter Q. Actually, Eq. lfTT]) . that defines our FW, is a superpo- 
sition of zero-order Bessel beams, and, due to this fact, the resulting field is expected to 
possess a transverse localization around p = 0. Each Bessel beam in superposition (JTTJ) 
is associated with a central spot with transverse size, or width, Ap n « 2A/k pn . On the 
basis of the expected expected convergence of series (jllj) . we can estimate the width of 
the transverse spot of the resulting beam as being 

2 4 2 4 

A P ~ T~ = / 2/2 n2 ' (16) 

which is the same value as that for the transverse spot of the Bessel beam with n = 
in superposition (JTTJ). Relation (fTtjj) can be useful: Once we have chosen the desired 
longitudinal intensity pattern, we can choose even the size of the transverse spot, and use 
relation Mb\) for evaluating the needed, corresponding value of parameter Q . 
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In spite of the fact that the transverse spot size happens to be approximately equal to 
that of a Bessel beam with k p = \Juj 2 /c 2 — Q 2 , it may happen that the decay of the field 
transverse intensity for p > Ap is much faster than that of an ordinary Bessel beam! This 
happens when the desired field intensity presents a longitudinal width Az much smaller 
than L, i.e., Az « L, as we will see below. 

An illustrative example: 
Let us consider the situation in which, within the interval < z < L, the desired longi- 
tudinal intensity pattern is given by a well-concentrated peak, represented by expression 
(HSJ), with A = 0.632 pm (that is, u = 2.98 x 10 15 Hz), Q = 0.98uj /c, L = 0.01m, and 
Az = L/500. 

Figures 4(a), 4(b) and 4(c) show the resulting FWs obtained by using, in superposition 
0, iV = 100, iV = 250 and N = 300, respectively. We can see in the first case that 
N = 100 is not enough for yielding a good result. On the other hand, the second and 
third cases, with iV = 250 and N = 300, seem to reproduce the desired pattern very 
well, with no apparent difference between the two cases. However, Fig. 5 shows that the 
residual intensity for the third case is smaller than for the second one, confirming the 
previous conclusions. 

Figure 6(a) shows the transverse intensity pattern of the peak (in the plane z = L/2) 
for the case with N = 250. We can see that the value of the transverse spot width agrees 
very well with our estimate ()16|). which furnishes, in this case, the value Ap « 1.22 pm. 
From Fig. 6(b) one can visually evaluate the residual intensity of the transverse pattern: 
One can observe that the transverse decay is strong and much faster than that presented 
by Bessel beams. This figure too confirms our previous conclusions. 
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Figure 4: Comparison of the desired longitudinal intensity pattern (solid line) with those 
of the resulting FWs (dotted line), when using: (a) N = 100; (b) N = 250; (c) N = 300. 




Figure 5: (a) Longitudinal residual intensity of the considered FW with iV = 250. (b) 
The same with iV = 300. 




Figure 6: (a) Transverse intensity pattern for the peak of the considered FW with N = 
250. (b) The transverse residual intensity for this case. 
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5 Increasing the control on the transverse shape by 
using higher-order Bessel beams 

We have shown in the previous Section how to get a very strong control over the 
longitudinal intensity pattern of a beam using suitable superposition of zero-order Bessel 
beams. 

As already mentioned, due to the fact that the resulting beam must obey the wave 
equation, we cannot get a total three-dimensional control of the wave pattern; but we 
have shown that we can have some control on the transverse behavior: More specifically, 
we can control the transverse spot size through the parameter Q, which defines the values 
of the transverse wave numbers of the Bessel beams entering superposition ([lip. 

In this Section we are going to argue that it is possible to increase even more our 
control of the transverse shape by using higher-order Bessel beams in our fundamental 
superposition (JTTJ). Despite the method presented in this Section is not yet demonstrated 
in a rigorous mathematical way, it can be understood and accepted on the basis of simple 
and intuitive arguments. The basic idea is obtaining the desired longitudinal intensity 
pattern, not along the axis p = 0, but on a cylindrical surface corresponding to p = p' > 0. 
This allows one to get interesting stationary field distributions, as static annular structures 
(tori), or cylindrical surfaces, of stationary light (or electromagnetic or acoustic field), and 
so on, with many possible applications.** To realize this, let us initially start with the same 
procedure in the previous Section; i.e., let us choose some desired longitudinal intensity 
pattern, within the interval < z < L, and calculate the coefficients A n by using Eq. fT0|) . 
Afterwards, let us replace the zero-order Bessel beams Jo(k pn p), in superposition (fTTj) . 
with higher-order Bessel beams, J M (A; pn p), to get 

N 

*(j>,z,t) = e -^ 0t e iQz Yl A n J,(k pnP )e l ^ nz , (17) 

n=-N 
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with A n = (1/L) Jq F(z)exp(—i2Trnz/L)dz, and k pn = — {Q + lnnjVf. 

In superposition (|T7|). the Bessel functions J p {k pn p\ with different values of n, reach 
their maximum values at p = p' n , where p' n is the first positive root of the equation 
(d J p (k pn p)/dp)\ p > n = 0. The values of p' n are located around the central value p' n=0 , at 
which the Bessel function J p {k pn= op) assumes its maximum value. We can intuitively 
expect that the desired longitudinal intensity pattern, initially constructed for p = 0, 
will approximately shift to p — p' n=Q . We have found such a conjecture to hold in 
all situations explicitly considered by us. By such a procedure, one can obtain very 
interesting stationary configurations of field intensity, as the mentioned "donuts" and 
cylindrical surfaces, and much more. 

In the following example, we show how to obtain, e.g., a cylindrical surface of sta- 
tionary light. To get it, within the interval < z < L, let us first select the longitudinal 
intensity pattern given by Eq.(|15p. with l\ = L/2 — Az and li = L/2 + Az, and with 
Az = L/300. Moreover, let us choose L = 0.05 m, Q = 0.998 u /c, and use N = 150. 

Then, after calculating the coefficients A n as before, 

A n = y[ L F(z)e-* 2 -? nz dz, 

we have recourse to superposition f|17jl . In this case, we choose p = 4. According to the 
previous discussion, one can expect the desired longitudinal intensity pattern to appear 
shifted to p' ~ 5.318/fc pn= o = 8.47/zm, where 5.318 is the value of k pn= op for which the 
Bessel function J 4 (fc pn=0 p) assumes its maximum value, with k pn=0 = \Juq — Q 2 . The 
figures below show the resulting intensity field. 

Figure 7(a) depicts the transverse intensity pattern for z = L/2. The transverse peak 
intensity is located at p = 7.75 pm, with a 8.5% difference w.r.t. the predicted value of 
8.47/im. In Fig. 7(b) the transverse section of the resulting beam for z = L/2 is shown. 

Figure 8 depicts the three-dimensional pattern of such a higher-order FW. In Fig. 8 (a) 
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the orthogonal projection of its 3D pattern is shown, which corresponds to nothing but 
a cylindrical surface of stationary light (or other fields). In Fig. 8(b) the same field is 
shown, but from a different point of view. 




Figure 7: (a) Transverse intensity pattern at z = L/2 of the considered, higher-order 
FW. (b) Transverse section of the resulting stationary field for z = L/2. 




Figure 8: (a) Orthogonal projection of the three-dimensional intensity pattern of the 
higher-order FW depicted in Figs. 7. (b) The same field but under a different perspective. 

We can see that the desired longitudinal intensity pattern has been approximately ob- 
tained, but, as wished, shifted from p = to p = 7.75 pm: and the resulting field resembles 
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a cylindrical surface of stationary light with radius 7.75 /zm and length 238 yum. Donut-like 
configurations of light (or sound) are also possible. 



6 Generation of Frozen Waves 

In the previous Sections, we have shown how suitable superpositions of Bessel beams 
of the same frequency can provide impressive results: Namely, can produce stationary 
wavefields with high transverse localization, and with an arbitrary longitudinal shape, 
within a chosen space interval < z < L; that is, Frozen Waves with a static envelope. 
As we already mentioned, such waves are rather interesting, not only from the theoretical 
point of view, but also because of their great variety of possible applications, ranging from 
ultrasonics to laser surgery, and from tumor destruction to optical tweezers. 
But how to produce our FWs? Regarding the generation of FWs, one has to recall that 
superpositions (fTTj). which define them, consist of sums of Bessel beams. Let us also 
recall that a Bessel beam, when generated by a finite aperture (as it must be, in any real 
situations), maintains its nondiffracting properties till a certain distance only (called its 
field-depth), given by 



where R is the aperture radius and 9 is the so-called axicon angle, related to the longitu- 
dinal wave number by the known expression cos 8 = c(5/uj. 

So, given an apparatus whatsoever capable of generating a single (truncated) Bessel 
beam, we can use an array of such apparata to generate a sum of them, with the appro- 
priate longitudinal wave numbers and amplitudes/phases [as required by Eq.()ll|)]. thus 
producing the desired FW. Here, it is worthwhile to notice preliminarily that we shall be 



Z 



R 




tan 9 
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able to generate the desired FW in the the range < z < L if all Bessel beams entering 
the superposition (jllj) are able to reach the distance L resisting the diffraction effects. 
We can guarantee this, for instance, if L < Z min , where Z min is the field-depth of the 
Bessel beam with the smallest longitudinal wave number (3 n= -N — Q — ZftN/L, that is, 
with the shortest depth of field. In such a way, once we have the values of L, ujq, Q, N, 
from Eq. (jl8|) and from the above considerations it results that the radius R of the finite 
aperture has to be 



R > L k 



^ --1 (19) 



The simplest apparatus capable of generating a Bessel beam is that adopted by Durnin 
et al. [7], which consists in an annular slit located at the focus of a convergent lens and 
illuminated by a cw laser. Then, an array of such annular rings with the appropriate radii 
and transfer functions, able to yield both the correct longitudinal wave numbers^ and the 
coefficients A n of the fundamental superposition (fTTjl . can generate the desired FW. 
In the next Section we shall just consider such a simple apparatus, even if, of course, 
other powerful tools, like the computer generated holograms, may be used to produce the 
FWs. 



6.1 A very simple apparatus for producing FWs 

Let us work out an example, by having recourse to an array of the very simple Durnin et 

al.'s experimental apparata. 

Since 1987, let us repeat, it has been used a simple experimental mean for creating a 

Bessel beam, consisting in an annular slit located at the focus of a convergent lens and 

illuminated by a cw laser. Let us call 5a the width of the annular slit, A the wavelength of 
^Once a value for Q has been chosen. 
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the laser, and / and R the focal length and the aperture radius of the lens, respectively. 
On illuminating the annular slit with a cw laser with frequency u , and provided that 
condition 5a <C \f / R is satisfied, the Durnin et al.'s apparatus creates, after the lens, a 
wavefield closely similar to a Bessel beam along a certain depth of field. Within such field- 
depth, z < Rj tan#, and to p << R, the generated Bessel beam can be approximately 
written 

4>(p,z,t)nAJ (k p p)e i ' 3z e iu!ot (20) 
with A a constant depending on the values of a, f, u and Sa, 

and 

P 2 = 4~ k l- ( 22 ) 

c 

Thus, as Durnin et al. suggested, we can see that the transverse and longitudinal wave 
numbers are determined by radius and focus of slit and lens, respectively. Once more, 
let us recall also that the wavefield has approximately a Bessel beam behavior (when 
p « R), in the range 0<z<Z^Rf/a that we have called the field-depth of the 
Bessel beam. Our FWs are to be obtained by suitable superpositions of Bessel beams. So 
we can experimentally produce the FWs by using several concentric annular slits (Fig.9), 
where each radius is chosen in order yield the correct longitudinal wave number, and where 
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the transfer function of each annular slit is chosen in order to furnish the coefficients A 
of Eq.Q which are needed for the desired longitudinal pattern to be obtained. 



Annular slits 




a T aN l r n 



5a i 



Lens 



Bessel beam 
superposition 



Sa«M 



Figure 9: A set of suitable, concentric annular slits, as a simple means for generating a 
Frozen Wave. 

Let us examine all this in more details. Suppose we have 2N + 1 concentric annular 
slits with their radii given by a n , with —N < n < N . Along a certain range, after the 
lens, one will have a wavefield given by the sum of the Bessel beams produced by each 
slit, namely* 



TV 

*(p,z,t) = e~ iu > ot J2 A n T n Jo(k P np)e if3nZ , (23) 

n=~N 



T n being the transfer function of the n-th annular slit (which regulates amplitude and 

phase of the emitted Bessel beam, and is a constant function for each slit); while the A n 

'■The same apparatus could also be used to generate higher order FWs, when the zero-order Bessel 
beams in superposition (|23|l are replaced with higher order Bessel functions. Experimentally, it can be 
performed by angular modulation of the slits. 
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are constants depending on the characteristics of the apparatus: Namely depending, in 
general, on the values of a, /, Uq and 5a. It is possible to obtain a simple expression for the 
A n by making some simplifying, rough considerations |Hj. The transverse and longitudinal 
wave numbers are given by 



Tv ryn 7~ 



"pn 



^0 <h 
C f 



(24) 



and 



Pi 



(25) 



On the other hand, we know from the present theory that for constructing the FWs 
quantity (3 is to be given by Eq.flZJ): 



On combining Eqs.(|71 (|2^j) . (|23jl . one gets 



^0 O-n 

of 



(26) 



and, solving with respect to a n , 



a n = f 



\ 



2tt 



1 - — [Q + -r n 



(27) 



25 



Equation (|27|) yields the radii of all the annular slits that provide the correct longitu- 
dinal wave numbers, needed for the generation of the FWs. We may notice that the radii 
of the annular slits do not depend on the specific desired longitudinal intensity-pattern, 
and that many different sets of values for the radii are possible on making different choices 
for the parameter Q. 

Notice that the procedure is not yet finished. Indeed, once the desired longitudinal 
pattern F(z) has been chosen, one has necessarily to meet in Eq.fjUJ) the coefficients 
A n given by Eq. (fT0| : and such coefficients have to be the coefficients of Eq. (f23| . For 
obtaining them, it is necessary that each annular slit be endowed by the appropriate 
transfer function, which regulates amplitude and phase of the Bessel beam emitted by 
that slit. By using Eqs. (|TU| ([TT |l . (p3]l . we get the transfer function T n of the n-th annular 
slit to be 




Finally, with the radius of each annular slit given by Eqs. (j27|) and the transfer func- 
tions of each slit given by Eqs. (}2*%j) . we do obtain a FW endowed with the desired longi- 
tudinal behaviour, inside the interval < z < L. Of course, one has to guarantee also 
that the distance L is smaller than the smallest field-depth of the Bessel beams entering 
superposition (|2*H|) . In other words, one must have also 



L < Z min « **J- (29) 
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where a max is the largest radius of the concentric annular slits. 



7 Conclusions 

In this work we have expounded the theory of Frozen Waves, and depicted some pos- 
sible experimental apparata to generate them. The present results can find applications 
in many fields: Just to make an example, in optical tweezers modelling, since we can 
construct stationary optical (but also acoustic, etc.) fields with a great variety of shapes; 
capable, e.g., of trapping particles or tiny objects at different locations.** These topics 
will be reported elsewhere. 
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